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Abstract 

We show that spontaneous breaking of global B — L symmetry responsible 
for small neutrino masses via the seesaw mechanism provides a unified picture 
of hybrid inflation and dark matter if the scale of B — L breaking is close to 
the GUT scale. The majoron which acquires a small mass in the milli-eV 
range due to Planck scale breaking of B — L is the dark matter candidate. 
The coupling of the majoron field to the neutrino induces a small violation 
of CPT and Lorentz invariance at the present epoch. We discuss some of the 
phenomenological and cosmological implications of the model. 
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I. INTRODUCTION 



While a non-zero neutrino mass is at the moment the only experimental evidence for 
physics beyond the standard model, as far as particle physics is concerned, in the domain of 
cosmology, several other oquestions remain unanswered within the standard model frame- 
work and also cry out for new physics. They are: (i) particle physics candidates for dark 
matter, (ii) a deeper understanding of the mechanism of inflation [1], (Hi) the nature of 
matter anti-matter asymmetry responsible for the observable universe [2] as well as (iv) 
an explanation of dark energy. There are many interesting and compelling models of new 
physics that can explain the different cosmological phenomena listed above individually. 
However, it is always much more desirable that a single model explain more than one phe- 
nomenon (ideally, of course, all of them). In this paper we discuss a model [3,4] originally 
designed to explain neutrino masses via the seesaw mechanism that seems to provide, after 
supersymmetrization, an explanation of both inflation and dark matter in a rather novel 
manner. The model has other desirable features such as gauge coupling unification, stable 
dark matter, stable proton due to R-parity conservation, etc. 

We use the supersymmetric singlet majoron [4] model described in Ref. [3] where the 
seesaw mechanism [5] for small neutrino masses is implemented by breaking a global B — L 
symmetry. The idea is to extend the standard model by the addition of one right handed 
neutrino per family and three standard model singlet superfields S, A and A. Of the new 
Higgs fields, S has B - L = 0, A has B - L = -2 and A has B - L = +2. The theoretical 
rationale for B — L symmetry is rooted in the present neutrino oscillation results which 
require that the seesaw scale be much lower than the Planck scale. It is therefore natural 
to think that it is protected by some symmetry. The simplest symmetry that does this 
is the B — L symmetry. B — L symmetry could be a global or a local symmetry. If we 
take it to be a global symmetry, its spontaneous breaking leads to a Nambu-Goldstone 
boson, the majoron. Since it is natural to expect that all global symmetries are broken 
by nonperturbative gravitational (or stringy) effects, we will parameterize these symmetry 
breaking effects by nonrenormalizable terms in the effective Lagrangian which are suppressed 
by the Planck scale M Pl . These Planck scale B — L breaking effects then give a tiny mass 
to the majoron [6]. 

In this paper we show that this model leads to: (i) a picture of hybrid inflation of the 
universe that links the neutrino mass (more precisly the seesaw or B — L scale of order 10 15 
GeV) to inflation with all the desired features and (ii) Planck scale B — L breaking effects 
that lead to a mass for the majoron in the milli-electron volt range and make it a candidate 
for dark matter of the Universe in the same way as the familiar axion, even though the 
parameters of the model are very different. It is, of course, interesting that this model ties 
the neutrino mass and the scale of inflation to the scale of gauge coupling unification. 

An interesting implication of the model is that it leads to a cosmologically induced 
Lorentz violation for neutrinos. We also find that the superpartners of the majoron must be 
in the MeV range in order to be compatible with the successes of big bang nucleosynthesis 
(BBN). 
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II. THE MODEL 



As already noted, the supersymmetric singlet majoron model [3] consists of the follow- 
ing superfields in addition to the well known quark, lepton, Higgs and gauge fields of the 
minimal supersymmetric standard model (MSSM): a right handed neutrino field N? (i is 
the generation index), new Higgs fields A, A which carry lepton number B — L = ±2 and 
a singlet field S, which is B — L neutral. We will show below that this model with an 
appropriate choice of the superpotential given below leads to F-term inflation as well as to 
the majoron as the dark matter candidate. We require the renormalizable part of the theory 
to be invariant under global U(1)b-l or an effective global symmetry which contains B — L 
as a subgroup 1 . 

The superpotential for the model can be written as a sum of three terms: 

W = Wmssm + Wo + Wl (1) 

where Wmssm is the familiar superpotential for the MSSM; Wo is the renormalizable part 
involving the new fields of the model and has the following form: 

W = h u LH u N + fNNA + XS(AA-v 2 BL ). (2) 

We assume that nonperturbative Planck scale effects induce nonrenormalizable terms in the 
potential and since we expect them to vanish in the limit of vanishing Newton's constant, 
they should be suppressed by powers of M P £. We arrange our theory such that the leading 
Planck scale induced term has the form 

Wi = -^{H u H d f{Af (3) 

At the nonrenormalizable level, there are many terms that one could write but we find this to 
provide a good description of physics of interest here and such a form could be guaranteed 
by additional discrete symmetries. For instance a symmetry under which S — > S; A — > 
iA; A -> -iA; H u -> H u ; H d -> iH d ; v c -> e^V; L -> e ilT / 4 L; e c -> e 5i?r / 4 e c ; d c -> 
— id c allows all the required terms for MSSM except the //-term. We include the /x-term in 
the superpotential since it is a lower dimensional term and it will break the symmetry softly. 

This superpotential leads to a potential involving the fields S, A and A with the following 
form: 



For instance one could have a local B — L model like a supersymmetric model based on the gauge 
group SU(2)l x U(l)i 3R x U(1)b-l and have local B — L symmetry be broken by a multiplet (S) 
with B — L > 6. The effective low energy theory in this case has a global U(l) symmetry which 
behaves like global B — L. One can construct a D-term inflation [11] model based on this model 
with a pair of A and S fields (and their conjugates), if we choose a form for the superpotential 
W = XS(AA - Ml) + AiXSS and a Fayet-Illiopoulos term for the gauged U(1) B -l- We do not 
elaborate on this model here. It has all the properties of the model presented here, with somewhat 
different parameters. 
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V(A,A,S) = A 2 |AA-<| 2 + A 2 |S| 2 (|A| 2 + |A| 2 ) + V 88 b\ (4) 

where Vssb stands for the supersymmetry breaking terms such as m 2 s \S\ 2 and |m 2 |A| 2 + 
m 2 |A| 2 etc. where these supersymmetry breaking mass parameters are all in the TeV range. 
We have ignored the Planck scale B — L breaking terms since they are small and not relevant 
to the discussion of inflation given in the next section. The minimum of this potential 
corresponds to (A) = (A) = vbl- 

Once the B — L symmetry is broken, this leads to the seesaw mechanism for small 

2 

neutrino masses. Typically, a neutrino mass has the form m v ~ where m v D is of the 
same order as a typical fermion mass of the standard model and Mr = fvsh with / — 1. If 
we assume that the atmospheric neutrino oscillation is linked to the third generation quarks 
(as would it be plausible in a quark- lepton symmetric theory), then the Dirac mass would 
be about 100 GeV and this would give vbl — 10 15 GeV, which is close to the value preferred 
by the inflation picture discussed below. 

Because the vevs of A and A break the global B — L symmetry, in the absence of explicit 
B — L breaking terms (denoted by Vpi anc k = 0), we have a massless particle in the theory, 
the majoron, given by = (x~X) I where we have parameterized A = -^(vbl+p)^ x ^ Vbl 
and A = -j=(vbl + p)e l ^l'" BL . The potential for is flat due to the shift symmetry under 
which — > + a. Once the Planck scale terms are turned on, the potential V(<f)) loses its 
flatness and can generate a rolling behaviour for the field as we see below. 



A. Inflation 

This model has the ability to generate inflation in early stages of the Universe. This 
comes about due to the interplay among the fields S and A and A. The potential involving 
them is given by Eq.(2), which has the form required in the hybrid inflation scenario. To see 
this note that for, S > vbl, the minimum of the potential corresponds to A = A = 0. We 
assume that A (f^) > (^)- The potential is then dominated by the term V = \ 2 v BL 
and the universe undergoes an inflationary phase. The S field keeps rolling towards the 
potential minimum and inflation ends when the field S reaches the value S ~ vbl- However 
it is necessary for the potential to have small tilt to drive the field toward its global minimum. 
This is the hybrid inflation picture [7]. 

One way to generate a slope along the inflationary trajectory (i.e. the S direction) is to 
include the one loop radiative correction to the tree level potential [8]. This arises because 
supersymmetry is broken during inflation and there will be a mass splitting between the 
components in the chiral multiplet. The effective one loop radiative correction is 

AV = J2 ^(- 1 ) 2Jl ( 2J ^ + V^s) ln (4jf ) (5) 

Where the sum extends over all the spin states J; with field dependent mass M { and A R is a 
renormalization scale. The F s = ^ ^ term splits the components of the chiral multiplets 
A and A into a pair of two real scalar fields with masses m| = ^-(S 2 ± 2v BL ) and a Dirac 
fermion with mass squared m 2 F = ^-S 2 . Since during the inflation A = A = 0, the effective 
potential reads 
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Ve ff =X\v% L + ^[(S* - 2v 2 BL ) 2 In (^^) + (S* + 2„| L )' In (^kl) - 25* In (**)]} 

(6) 

For 5 much larger than the (B — L) symmetry breaking scale the above potential becomes 

I^AkO + ^ln^} (7) 

The log term will provide the driving force for the field S to roll down the potential. In this 
case the slow roll parameters are given by 

1 _ _P 1_ 

lvl ~ 2N(S)^ e ~ 32tt 2 N(S) ( ) 

where N is the number of e-foldings and is given by 

Density perturbations in about N eo = 60 e-foldings before the end of inflation are estimated 

as 



*p ^ i6 J ^(^r (io) 

P V 3 K M P / 1 ; 

Fluctuations with amplitude ~ 10~ 5 can obtained if vbl — 6 • 10 15 GeV. 

Another way to generate a slope to the potential is to add soft supersymmetry breaking 
mass ms to the inflaton field [9]. In this case the slow roll parameters are given by 

A M Pi m 2 .m 2 S 2 , 

^8^» £ (11) 
The requirement that rj <^ 1 (say 0.01) gives 

A ^ 10 3/ 2( ^ )( ^) (12) 
vbl vbl 

In order for the soft supersymmetry breaking to dominate over the one loop radiative correc- 
tion the parameter A must be smaller than 3 • 10~ 5 . The value of S after suffering N m ~ 60 
e-foldings between the horizon exit and the the end of inflation is given by 

5*60 = Vbl cxp (6O77) (13) 

The A field plays the role of the inflaton field, which then oscillates and reheats the Universe 
via the production of neutrinos due to NNA coupling. The density fluctuations lead to 
strong constraints on the parameters of the model and are given by : 
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^AVM^) 3 (^) 2 exp(-6 0?7 ); (14) 
p Mp£ ms 

For vbl — 2 • 10 14 GeV , m s ~ 2 TeV and A ~ 1CT 4 - 5 we get 5p/p ~ 1CT 5 . It is interesting 
that the value for vbl required for understanding neutrino masses also gives the correct 
order of magnitude for density fluctuations generated by inflation. 

Finally, we discuss the question of reheating after inflation ends. The reheating can 
be assumed to be caused by the decay of the field S when it starts oscillating around the 
minimum. The reheating temperature is given by the approximate formula T R ~ ^/MpiTs 
where is the decay width of the inflaton. To obtain the decay width of S, we have to 
isolate the decay modes. We expect that the A coupling in the renormalizable part of the 
superpotential will give the dominant contribution to the decay width. 

To see this in detail, note that if both A and A have equal vev 2 , the combination 
ip = -^(ipA— V'a) is the fermionic partner of the majoron field (denoted by 0) and is therefore 
light whereas ^(V'A + V'a) is the fermionic partner of the superheavy supermultiplet and 
therefore has same mass as S field. However one can see from the superpotential that S has 
a coupling to ipip and can therefore decay to these states. The decay width is then given by 
Ts ~ ^ m s — 10~ 14 ' 5 i>.bl- Note that the inflaton gets its dominant mass from the B — L 
breaking vev vbl and a smaller mass from supersymmetry breaking. Therefore M s — vbl- 
This leads to a reheating temperature of T R ~ 10 9 GeV. In this case a gravitino with mass 
m 3/ / 2 — lOOGeV does not spoil the success of big bang nucleosynthesis. 

It is also worth pointing out that the SUGRA embedding of this model os free of the 77 
problem that is generic to F-term SUGRA models [11]. 



B. Nonrenormalizable terms, majoron mass and majoron as the ultralight dark 

matter: 

In the presence of the nonrenormalizable terms, the majoron picks up a nonzero mass. If 
we choose the leading order non-renormalizable term to be y£lun *A y - ' which can be done 

by a judicious choice of discrete symmetries as discussed above), then we get the majoron 
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mass rr£ ~ — (KT 15 eV) 2 . 

In order to discuss how massive majorons constitute the dark matter of the universe, 
let us discuss the evolution of the majoron field as the universe evolves. The majoron 
potential, which arises solely from the B — L breaking terms can be written using Eq. (3) 
in the form: 

V{4>) = AiA 4 (l + cos^) , Ai~l (15) 

vbl 

where A]A 4 = m 2 ^v BL ~ (1.7 x 1CT 3 GeV) 4 . We have added a constant term, akin to a 
cosmo logical constant, so that after the amplitude for the <fi field oscillation damps to its 



2 these vevs will differ by a small amount proportional to the supersymmetry breaking which can 
therefore be neglected for large values of vbl that we are interested in. 
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minimum at = itvbl the value of the potential is zero. The field and the radiation 
energy density p satisfy the coupled scalar field - Einstein-Friedman equations 



<t> + 3H<j>=-V'(<f>) (16) 

and 

In the early universe, clearly the V term is negligible. The Hubble parameter H is then 
dominated either by the kinetic energy of or the relativistic energy density in radiation. 
Solving Eq. 16, one finds that oc R~ 3 where R is the scale factor. So regardless of whatever 
initial the field starts out with, it completely damps down to a very small value as the 
universe expands. Thus the value of freezes over for most of the early universe at some 
random value. When ~ 3H, the potential term will dominate Eq. (16) and the field will 
oscillate around its minimum value of = ttvbl- These oscillations, which do not damp, 
behave like matter [12] and contribute to the dark matter density. Below, we calculate this 
contribution and find that for the range of parameters of interest, it gives the right order of 
magnitude for Vt matter . 

As noted already, the value of the field remains frozen at its initial value (taken to be 

3T 2 

vbl) until the epoch when ~ 3H ~ pj- at which time the field starts to oscillate. The 

temperature at which it does this is given by Tj ~ ^Mpgm^jZ ~ 1CT 3 GeV. As it oscillates, 
it is easy to show that it acts like a pressureless gas and thus can be treated as an ensemble 
of nonrelativistic particles [12]. The contribution of the majoron to the energy density now 
is given by: 

The ratio of the majoron (dark matter) to the radiation energy density at their equipartition 
temperature T E q ~ 1 eV is given by 

* (19) 

PR p R lEQJ-i 

where T ~ 2.4 • 10~ 4 eV and p^ ~ 2 • 10~ 15 eV 4 are the temperature and the energy density 
of the photons today. For Ax A 4 ~ 1 MeV 4 and Tj ~ 100 MeV one obtains ^(T EQ ) ~ 1. 
Hence the universe remains radiation dominated until the temperature drops below ~ 1 eV. 

Thus we see that majoron oscillations can act as a nonrelativistic dark matter with the 
right order of magnitude for fl ma tter- We therefore conclude that the same range of parame- 
ters that gave the neutrino mass values in the right order and also the density fluctuations 
required by CMB measurements also seems enable the majoron to play the role of dark 
matter. 

To make these ideas concrete, we have solved numerically the coupled equations (16) 
and (17) in the regime prior to and after the majoron field begins oscillating. That is the 
region in which the Hubble parameter, which includes the sum of radiation, matter and the 
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majoron field energies, in the second term in the right hand side of Eq. (16), has decreased 
to the point that the right hand side has become appreciable in comparision with that term. 
For natural choices of the parameters discussed above (vbl and v wk 10~ 4 Mp£ and 100 GeV 
respectively) the coefficient A 4 in the Eq. (15) becomes (10~ 4 GeV) 4 . In Figures 1 and 2, 
we plot the solution for the majoron field for two initial values: 4>/vbl = 0.17T and ir/3. 
We note that the amplitude of the final field oscillation is less in the second case. For still 
larger initial values of 0, the amplitude falls off more sharply. On the other hand, both for 
smaller initial values and over a wide range of values, the final result for the majoron 
oscillation amplitudes are the same. 

In Figures 3 and 4 we plot the resulting ratio p^/ p 1 for the same choice of parameters as 
in Figures 1 and 2 respectively, over a broad region in which the oscillatory frequency is not 
too large compared to the expansion rate so that numerical computations with reasonable 
number of steps can be carried out and are reliable. One sees that this ratio increases 
linearly with decreasing temperature, as is appropriate for presureless matter. One can find 
today's value for p^j p 1 simply by multiplying the values read-off the figure by T/T . These 
numerical computations confirm that our parameter choice leads to the correct order of 
magnitude for today's ratio (that is, oscillating field giving Qdm ~ 0.3). 

Finally we note that the frequency of the majoron field oscillation is given by uo = 
(XiA A /v BL ) 1 ^ 2 ~ 0.2 sec -1 . The oscillations begin around the epoch of nucleosynthesis. The 
energy in these oscillations around this epoch is very small compared to that in p 7 so that it 
does not affect the BBN considerations. Furthermore it is worth noting that the frequency 
uj is independent of the value of vbl- 



III. OTHER CONSEQUENCES OF THE MODEL 



Some consequences of the superlight majoron for cosmology has recently been discussed 
in [13]. We consider other implications that particularly relate to the model in this paper 
in this section. 



A. Cosmologically induced Lorentz violation for neutrinos 

One important implication of our proposal is that it leads to a cosmologically induced 
violation of Lorentz invariance. This comes about because in our model the majoron field has 
a derivative coupling of the form -^-uj^ud^cp. In the late universe when ^ 0, this leads to 

an effective Lorentz violating term of the form -^y^v in the effective low energy Lagrangian. 
This effect is Lorentz violating and will manifest itself in the neutrino oscillation process 
[14]. The maximum value of the parameter in our model is m$ ~ 10~ 23 — 10~ 24 GeV. 
However since the Lorentz violating term is family universal, it will most likely manifest irself 
in a transition from v a to Vp. So the only way to detect this will be to measure P[y a — > 
up) — P(P/3 — > u a ). In the Lorentz invariant case, such effects are generally suppressed by the 
mass of the neutrino [15] and are therefore likely to be very small. The detailed experimental 
implications of such Lorentz violating effects are currently under investigation. 
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B. BBN constraints on majoron superpartners 



A second phanomenological implication is that if we ignore the supersymmetry breaking 
effects, then majoron belongs to a supermultiplet together with a scalar partner (to be 
called smajoron) and a fermion, majorino, both of which are massless. In the presence 
of supersymmetry breaking effects, the smajoron (a) and majorino (-0) P^k U P mass in 
the MeV to GeV range. The precise values of these masses are however constrained by 
cosmological consideration for the model to be viable. We discuss them in this section. 

This question was discussed in [3] for the case where the B — L breaking scale is in the 
TeV range. It was found that in that case the masses of a and tp can be in the TeV range. 
This holds as long as vbl < 10 8 GeV. In our case however vbl — 10 15 GeV. We will therefore 
end up with the smajoron (er) and majorino (tp) masses in the 10 MeV or lower range. 

The primary goal is to make sure that the new particles do not affect big bang nucle- 
osynthesis (BBN). Since the couplings of the 0, a and ip to standard model particles are 
all suppressed by vbl, the lifetimes of a and -0 are much longer than one second, the BBN 
epoch. Therefore if they are heavy and their abundances at the BBN epoch are not sup- 
pressed, then they will have adverse effect on nucleosynthesis. We therefore have to calculate 
their abundance at the BBN epoch. 

Due to the suppressed couplings of 0, a and -0, they will decouple in the very early 
stage of the universe. To calculate the decoupling temperature, we note that the typical 
annihilation rates are for 0, a and -0 are given by R ~ Using the decoupling condition 

V BL 

R{T D ) < H(T D ), we get 

t d - 9l /2 (T D )^ (20) 

which leads to To — 10~ 1 vbl which is of order 10 14 GeV. After decoupling the 0,(7 and 
■0 densities simply dilute due to the expansion of the Universe. Their contribution to the 
energy density of the universe at the BBN epoch is given by 

P<t>,a,i(> ^ n 4>,v,'4> m <i>,<T,i> ^ 9*{Tbbn) m<j,,u^ 
Pt n y Tbbn 9*{Td) Tbbn 

Success of BBN requires that should be much less than one. Since 9 * ( - T , I L B ^ ~ 100, we 

^ Pi 9*(Td) 

therefore must have <C 100 MeV. 

We further note that since these particles 0, a, -0 are singlets under the standard model, 
they are not in conflict with any known low energy observations despite their small mass. 



C. Implications for leptogenesis 

Finally, we wish to comment that in this model one can employ the mechanism of lep- 
togenesis [16] to understand the origin of matter in the universe. One can have a right 
handed neutrino [17] at the intermediate scale range of 10 9 GeV range (or even a pair of 
them nearly degenerate [18]), whose decay would produce a lepton asymmetry, which via the 
sphaleron interactions would get converted to baryons. The only new aspect of our model 
is the presence of the majoron at very low energies. In principle its interactions can erase 
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the baryon asymmetry since its interactions violate lepton number. However in our case 
since the scale vbl is very high, as noted already, the majorons decouple around 10 14 GeV 
and are therefore "impotent" as far as their effect on lepton asymmetry is concerned. If for 
instance the vbl scale was in the range below 10 9 GeV, we would have no leftover lepton 
asymmetry at the weak scale to be converted to baryons. It is therefore interesting that the 
high scale of vbl is required from various considerations. 

In conclusion, we have presented a simple model for neutrinos using the super symmetric 
extension of the singlet majoron model that provides a unified framework for understanding 
inflation and dark matter for the same set of parameters required by neutrino masses. We 
find that the vbl scale is constrained from various considerations to be around the conven- 
tional grand unification scale ~ 10 15 GeV. We have checked our results using a numerical 
solution of the evolution equation for the majoron field. The high vbl scale makes the 
majoron and its superpartners highly invisible in collider experiments. 

The works of R. N. M. and S. N. are supported by the National Science Foundation 
Grant No. PHY-0099544. The works of D. K. and V. L. T. are supported by NASA. R. N. 
M. would like to thank Z. Chacko for discussions. 
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FIG. 1. The evolution of the Majoron field (in units of Mpi) as a function of temper- 
ature for the initial value of the field <^>(0) = 10~ 5 . The field settles down to the value 
4> = kvbl', vbl = W~ 4 Mpi, driven by the first minumum of the potential. The Log is with 
respect to base 10 and e -4 in the ordinate stands for 10~ 4 . 
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FIG. 2. The evolution of the Majoron field (in units of Mpi) as a function of temperature for 
the initial value of the field (f)(0) = 10~ 4 . In this case also, the field settles down to the same value 
(f) = ttvbl', v bl = 10 Mpi, driven by the first minumum of the potential. The notation is same 
as in Figure 1. Note that the amplitude of oscillation is smaller than that in Fig. 1 by about a 
factor of half. 
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FIG. 3. The ratio of energy density in the field (f>, p§ = <ft 2 /2 + V(</>) over the radiation energy 
density p ra( i as a function of the temperature of the Universe for initial value of 0(0) = 10~ 5 Mp;. 
After the field begins oscillating it behaves just like matter, as it can be seen from the above ratio 
which increases like T. 
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FIG. 4. The ratio of energy density in the field (f>, ps = (j) 2 /2 + V(4>) over the radiation energy 
density p ra( i as a function of the temperature of the Universe for initial value of (p(0) = 10 _4 Mp;. 
After the field begins oscillating it behaves just like matter, as it can be seen from the above ratio 
which increases like T. 
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